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ABSTRACT 

KA1858+4850 is a narrow-line Seyfert 1 galaxy at redshift 0.078 and is among the brightest active 
galaxies monitored by the Kepler mission. We have carried out a reverberation mapping campaign 
designed to measure the broad-line region size and estimate the mass of the black hole in this galaxy. 

We obtained 74 epochs of spectroscopic data using the Kast Spectrograph at the Lick 3-m tele- 
scope from February to November of 2012, and obtained complementary V-band images from five 
other ground-based telescopes. We measured the H/3 light curve lag with respect to the V-band 
continuum light curve using both cross-correlation techniques (CCF) and continuum light curve vari- 
ability modeling with the JAVELIN method, and found rest-frame lags of tccf = 13.53^2 32 days and 
UjAVELIN = 13.151^qq days. The H/3 root-mean-square line profile has a width of <Tii ne = 770 ± 49 km 
s _1 . Combining these two results and assuming a virial scale factor of / = 5.13, we obtained a virial 
estimate of Mbh = 8. 06+]/ 72 x 10 6 M 0 for the mass of the central black hole and an Eddington ratio 
of L/I/Edd ~ 0.2. We also obtained consistent but slightly shorter emission- line lags with respect to 
the Kepler light curve. Thanks to the Kepler mission, the light curve of KA1858+4850 has among 
the highest cadences and signal-to- noise ratios ever measured for an active galactic nucleus; thus, our 
black hole mass measurement will serve as a reference point for relations between black hole mass and 
continuum variability characteristics in active galactic nuclei. 
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1. INTRODUCTION 

The NASA Kepler Mission, designed to search for exo- 
planets, continuously monitored the brightness of more 
than 100, 000 stars in a 115 square-degree field for about 
four years (Borucki et al. 2010). Situated within the Ke- 
pler field are several active galactic nuclei (AGNs) that 
also exhibit optical flux variations. Kepler's monitoring 
capabilities enable measurements of AGN optical light 
curves over long temporal baselines with unprecedented 
cadence and precision, providing the basis for extremely 
detailed AGN variability studies. 

Observations have revealed correlations be- 
tween AGN variability amplitude and redshift 
(Cristiani et al. 1990; Giallongo et al. 1991; Hook et al. 
1994; Cid Fernandes et al. 1996; Vanden Berk et al. 
2004), variability amplitude and black hole mass 
(Wold et al. 2007; Wilhite et al. 2008; Bauer et al. 

2009) , and anticorrelations between variability ampli- 
tude and luminosity (Cristiani et al. 1990; Hook et al. 
1994; Cid Fernandes et al. 1996; Giveon et al. 1999; 
Vanden Berk et al. 2004; Webb & Malkan 2000). Fur- 
thermore, analyses of continuum light curves have 
revealed the presence of characteristic variability 
timescales which have been found to vary with black 
hole mass (Collier & Peterson 2001; MacLeod et al. 

2010) . The Kepler high-resolution light curves have a 
cadence of 30 minutes, and are the only datasets to date 
that have been able to probe optical AGN variability 
down to such short time scales. Optical fluctuation 
power spectral density functions for several Kepler 
AGNs have already been published (Mushotzky et al. 

2011) , and they have shown much steeper slopes than 
those seen in the X-rays. Kepler's light curves provide 
new high signal-to-noise ratio (SNR) data which will 
test and better constrain these previously established 
correlations and further shed light on AGN variability 
characteristics. 

Independent measurements of black hole mass are re- 
quired to search for connections between AGN vari- 
ability characteristics and black hole mass. To this 
end, we present the results of a nine-month monitoring 
campaign for the narrow- line Seyfert 1 (NLS1) galaxy 
1RXSJ185800. 9+485020, also known as KA1858+4850, 
which has redshift z = 0.078 and a Galactic extinc- 
tion of Ay = 0.15 mag (Schlafly & Finkbeiner 2011). 
This object was identified as an X-ray source in the 
ROSAT All-Sky Bright Source Catalogue (Voges et al. 
1999). Prior to 2012, there was no published spec- 
trum of KA1858+4850 in the literature, and an obser- 
vation from Lick Observatory identified it as a Seyfert 
1 galaxy (Edelson & Malkan 2012). The initial portion 
of KAl858+4850’s Kepler light curve from quarters Q6 
and Q7 was published by Mushotzky et al. (2011) and 
showed strong optical variability, qualifying it as a prime 
candidate for reverberation mapping. 

The technique of reverberation mapping relies on the 
assumption that variability in the AGN continuum is 
echoed by emission lines originating from the surround- 
ing broad-line region (BLR; Blandford & McKee 1982). 
Ionizing photons from the AGN central engine travel to 
the BLR gas in a time r that is a function of the BLR 
radius. Changes in the ionizing photon flux incident on 
BLR clouds cause fluctuations in the emission-line flux. 


This means that the emission-line light curve will appear 
as a lagged version of the continuum light curve, and 
the lag time, combined with the speed of light, can give 
an estimate of the BLR radius. Additionally, the line- 
emitting gas orbits the central black hole at very high 
velocities, which causes Doppler broadening of the emit- 
ted spectral lines. The width of the broad emission line 
gives the velocity dispersion of the BLR gas, which, com- 
bined with the BLR radius, can yield a virial estimate of 
the central black hole mass. 

Kepler light curves covering over two years of moni- 
toring are now publicly available for KA1858+4850, of 
which three consecutive quarters (Q13, Q14, and Q15) 
directly coincide with the time of our ground-based mon- 
itoring campaign. We therefore performed our analysis 
using both F-band and Kepler observations. 

We employed the Lick Observatory 3-m Shane tele- 
scope with the Kast Spectrograph and five other ground- 
based telescopes to spectroscopically and photometri- 
cally monitor KA1858+4850 from February to Novem- 
ber of 2012. We describe our imaging observations and 
data reductions in §2 and §3; spectroscopic observations, 
reductions, and measurements are described in §4 and 
§5; §6 outlines the steps in measuring emission-line light 
curve lags; our estimates of the black hole mass (Mbh) 
and Eddington ratio are discussed in §7 and §8; and §9 
summarizes our results. 

2. IMAGING OBSERVATIONS 

Reverberation mapping requires a continuum light 
curve with high sampling cadence and SNR. To achieve 
this, we obtained F-band images from ground-based 
telescopes and used aperture photometry to construct 
a light curve for KA1858+4850 that has nearly nightly 
sampling for a span of 290 days. For several reasons, we 
chose to use the F-band light curve rather than the Ke- 
pler light curve for reverberation measurements. First, 
we wanted to monitor the AGN’s variability in real time, 
and since Kepler data are uploaded only periodically, this 
was possible only with ground-based monitoring. Addi- 
tionally, the Kepler passband, at 4000—8650 A, includes 
the strong Ha emission line, which can contribute signif- 
icantly to the photometric fluxes and introduce a strong 
lag signal to what should ideally be a pure continuum 
light curve. Furthermore, Kepler light curves exhibit 
severe mismatches between the flux scales for different 
quarterly observation sets, as can be seen in light curves 
shown by Revalski et al. (2014). We avoided these is- 
sues by constructing the continuum light curve with pho- 
tometric data from five ground-based telescopes, whose 
properties are described in the following sections. 

2.1. West Mountain Observatory 

The Brigham Young University West Mountain Ob- 
servatory (WMO) uses a 0.9-m telescope that employs 
a FLI PL3041UV detector with a 2018 x 2018 field of 
view. The CCD has 15 pm pixels and a scale of 0'.'61 
pixel -1 . KA1858+4850 was observed at WMO with ex- 
posure times of 200 s, 240 s, 250 s, or 300 s. WMO data 
covered the period from March to November of 2012 with 
images from 124 nights, and had a median seeing of 3'/ 2. 
Figure 1 shows a portion of the WMO held of view cen- 
tered on KA1858+4850. 
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2.2. KAIT 

The Katzman Automatic Imaging Telescope (KAIT) 
at Lick Observatory is a 0.76-m robotic telescope with an 
Apogee AP7 CCD, which has 24 ^m pixels in a 500 x 500 
array and a scale of 0'/8 pixel -1 (Filippenko et al. 2001). 
KA1858+4850 was observed with KAIT using 300 s ex- 
posures with the exception of six nights, for which expo- 
sure times of 60 s, 180 s, or 240 s were used. The median 
seeing for the KAIT exposures was 3"2, and the observ- 
ing period at KAIT spanned February to September of 
2012 with data from 109 nights. 

2.3. Faulkes Telescope North 

The Faulkes Telescope North (FTN), operated by the 
Las Cumbres Observatory Global Telescope Network, is 
a 2-m telescope located at the Haleakala Observatory in 
Hawaii. We used the Spectral camera with a Fairchild 
Imaging CCD486 detector, which has a 10f5 x 10^5 field 
of view (Brown et al. 2013). The CCD has 15 /im pixels 
in a 4000 x 4000 array and has a scale of 0 , . , 152 pixel -1 . 
The images were obtained using 2x2 binning for the 
readout. KA1858+4850 was observed at FTN with 120 
s exposures from February to March 2012. The exposure 
time was increased to 180 s in April 2012, then to 240 
s in May for the remainder of the program ending in 
November. We obtained 65 epochs of data from FTN, 
with median seeing of l'/6. 

2.4. The Nickel Telescope 

The 1-m Nickel telescope at Lick Observatory employs 
a Loral 2048 x 2048 CCD with a 613 x 613 field of view 
and a scale of 0'/184 pixel -1 . The images were obtained 
using 2x2 binning for the readout. KA1858+4850 was 
observed on the Nickel adopting 300 s exposures with 
the exception of three nights, for which 150 s, 250 s, and 
600 s exposures were used. We obtained 47 epochs of 
data from the Nickel between February and November of 
2012, and the median seeing was 2'/4. 

2.5. Mount Laguna Observatory 

The Mount Laguna Observatory (MLO) 1-m tele- 
scope uses a Fairchild CCD that has 15 /m i pixels in 
a 2048 x 2048 array, and has a scale of 0'/41 pixel -1 . 
KA1858+4850 was observed at MLO with 300 s expo- 
sures. The median seeing for the MLO exposures was 
3'.'0. Between February and November of 2012, we ob- 
tained 27 epochs of data from MLO. 

3. PHOTOMETRIC REDUCTIONS AND MEASUREMENTS 
3.1. V-Band Data 

Photometric data reduction included overscan cor- 
rection, trimming, bias subtraction, and flat fielding 
for all images. We used the Astrometry.net software 
(Lang et al. 2010) to register celestial coordinates onto 
images from WMO, KAIT, Nickel, and MLO. This step 
was omitted for FTN data, which already contained ce- 
lestial coordinates in the image headers. After cleaning 
all images of cosmic rays using the L. A. Cosmic algorithm 
(van Dokkum 2001), we performed aperture photometry 
in IDL using an aperture radius of 3" and sky annulus 
radii of 10"— 20", and obtained instrumental magnitudes 
for KA1858+4850 and seven comparison stars (marked 



Fig. 1. — A subset of a coadded frame created from WMO im- 
ages showing KA1858+4850 (boxed) and its seven comparison stars 
(circled). 


in Figure 1) for each image. The comparison stars were 
chosen to have similar or slightly brighter F-band mag- 
nitudes compared to KA1858+4850. For nights where 
multiple exposures were taken at the same telescope, the 
magnitude measurements for each object were averaged 
into a single value. Since KA1858+4850 is almost indis- 
tinguishable from a point source at ground-based reso- 
lution, we did not attempt to remove host-galaxy light 
from the AGN photometry. 

We used the comparison stars as constant-flux ref- 
erences and obtained a separate AGN light curve for 
each telescope. However, the uncertainties from aper- 
ture photometry photon counting errors underestimate 
the true photometric error budget. Additional sources of 
error include inconsistencies in flat-field corrections and 
poor comparison-star magnitude measurements owing to 
blemishes on the detector. We measured the magnitude 
of these additional errors by calculating the excess vari- 
ance, defined as 

a l = T? - ( 1 ) 

i=l 

in the scaled comparison-star light curves. Here, N is the 
number of measurements in the sample, /i is the mean 
magnitude, and Xi and cr,; are the individual measure- 
ments and their associated uncertainties, respectively. 
The cr,; values range from 0.004 mag to 0.048 mag, and 
the median and standard deviation of the uncertainties 
are 0.009 mag and 0.007 mag, respectively. We found 
the mean scatter of all seven comparison stars to be 
cr x « 0.001 mag, and added this in quadrature to the 
uncertainties from aperture photometry to produce the 
final AGN light curve for each telescope. 

To combine the light curves from different telescopes, 
we scaled each light curve so that the mean comparison- 
star magnitudes for each telescope matched those from 
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TABLE 1 

Photometric Comparison Stars for KA1858+4850 


Star 

a 

( h:m:s ) 

<5 

(°:':") 

V 

(mag) 

i 

18:58:04.03 

48:51:15.53 

16.612 ± 0.026 

2 

18:58:06.24 

48:51:19.82 

17.121 ± 0.037 

3 

18:58:09.54 

48:50:20.69 

17.256 ± 0.031 

4 

18:58:11.03 

48:49:46.16 

15.449 ± 0.028 

5 

18:58:09.97 

48:48:54.06 

15.195 ± 0.029 

6 

18:58:06.09 

48:49:16.65 

15.843 ± 0.029 

7 

18:57:58.01 

48:51:49.72 

16.724 ± 0.027 


Note. — Coordinates are J2000 and are based on an astromet- 
ric solution obtained by the astrometry.net software (Lang et al. 
2010). The quoted uncertainties are calculated as the standard 
deviation of 18 measurements from photometric nights at WMO. 

WMO, the telescope with the highest SNR and cadence 
and longest temporal coverage. However, each telescope 
has a different wavelength-dependent throughput, which 
can cause systematic offsets between light curves from 
different telescopes since the AGN is likely bluer than the 
average comparison-star color. We tested for these offsets 
by calculating the differences between AGN magnitude 
measurements taken on the same night but at different 
sites, and found the offsets to be on the order of 0.01 
mag. We applied these calculated shifts to the FTN, 
KAIT, MLO, and Nickel light curves and brought them 
into agreement with WMO to produce the combined light 
curve. 

Finally, we used Landolt (1992) standard stars ob- 
served at WMO to calibrate the zero point of the mag- 
nitude scale and produce the final light curve. We 
used WMO images from 18 nights, on which the ob- 
servers deemed conditions photometric, to calibrate the 
comparison-star magnitudes. We did not attempt to 
compute color dependence in the Landolt calibrations. 

Because truly photometric conditions are rare and 
difficult to confirm, each night gave slightly different 
comparison-star magnitudes. We took the weighted 
mean magnitude and standard deviation over 18 nights 
to be the magnitude and uncertainty for each star. The F 
magnitudes of the comparison stars are listed in Table 1. 

Figure 2 plots the final F-band light curve for 
KA1858+4850. The vertical length at each epoch in- 
dicates the photometric uncertainties, and the data are 
listed in Table 2. We averaged photometric measure- 
ments taken within 12 hours of each other to produce a 
condensed light curve that was used for subsequent lag 
analyses. 

The steps from performing aperture photometry to ob- 
taining a multiple-telescope light curve were carried out 
using an automated pipeline. Mapping WCS coordinates 
onto the images allowed for automatic detection of the 
AGN and comparison-star locations for aperture pho- 
tometry. The automated nature of this process enables 
the pipeline to process a large number of images at once, 
and to rapidly produce and update the AGN light curve 
as new images are acquired. 

3.2. Kepler Data 

We also obtained Kepler Simple Aperture Photometry 
(SAP) fluxes for KA1858+4850 from the MAST archive 
for Q13, Q14, and Q15, corresponding to 2012 March 
through November. Data from Q12 are missing from 
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Fig. 2. — KA1858+4850 E-band light curve. The vertical length 
at each epoch indicates the photometric uncertainties. 


TABLE 2 

Photometry measurements for KA1858+4850 


UT Date 

Telescope 

HJD - 2450000 

V (mag) 

2012-01-31 

N 

5957.754 

17.129 ± 0.008 

2012-02-11 

M 

5969.014 

17.012 ± 0.015 

2012-02-17 

N 

5974.663 

17.055 ± 0.015 

2012-02-20 

F 

5978.153 

17.011 ± 0.014 

2012-02-22 

F 

5980.166 

17.037 ± 0.020 


Note. — The telescopes are listed as follows: N = Nickel, M = 
MLO, F = FTN, K = KAIT, W = WMO. (Full table available in 
online version.) 


the archive because, during this time, the source fell on 
Module 3 of the Kepler telescope, which failed early on 
in the mission. 

The Kepler light curves are mismatched between indi- 
vidual quarters, so we used our F-band light curve as a 
reference to scale each quarter’s light curve individually. 
We applied a different multiplicative scale factor and ad- 
ditive shift to each quarterly Kepler light curve to bring 
it into agreement with ground-based observations. The 
multiplicative factors account for the difference in trans- 
mission between the Kepler and F-band filters, and the 
additive constants account for the changes in AGN-to- 
host galaxy flux ratio between each quarter caused by 
using different quarterly extraction apertures to obtain 
SAP fluxes. 

For each epoch in the condensed F-band light curve, 
we averaged together all Kepler flux measurements taken 
within six hours of the F-band measurement to compose 
condensed Kepler light curves. Then for each quarter, 
we fitted the contemporaneous Kepler and F-band flux 
measurements to the equation 

fv = m* /Kepler + b , (2) 

where m gives the multiplicative scale factor and b gives 
the additive shift. We fitted the data using MPFITEXY 
to account for measurement errors in both F-band and 
Kepler data. Figure 3 shows the results of applying a 
scale factor (top panel) and a scale factor plus a shift 
(middle panel) to the Kepler light curves. 

Even with a multiplicative scale factor and an addi- 
tive shift, however, there are still visible discrepancies 
between the two sets of data. Specifically, each Kepler 
quarterly light curve tilts downward with time compared 
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Fig. 3. — F-band light curve (black) overplotted with Kepler 
Q13, Q14, and Q15 light curves in red, blue, and green (respec- 
tively), scaled by three different methods. Top: multiplicative 
factor only; middle: multiplicative factor plus additive constant; 
bottom: multiplicative factor plus additive constant fitted to a V- 
band light curve with an additional linear trend. Error bars are 
not plotted for the Kepler points. 


to the F-band light curve. This is caused by the con- 
stant change in Kepler pointing with respect to the Ke- 
pler field as the telescope orbited the Sun, which, in turn, 
causes differential velocity aberration (DVA) and results 
in a trend that is superimposed on the light curve within 
the period of each quarter (Still & Barclay 2012). To 
account for this effect, we applied an empirical secular 
linear trend to the F -band light curve by adding a time- 
dependent flux to the data. The Kepler light curves were 
then fitted to the adjusted F-band light curve with scale 
factors and shift constants. Finally, the empirical trend 
was removed from both F-band and Kepler light curves 
by subtracting the same time-dependent fluxes as before. 
The resulting scaled Kepler light curves are shown in the 
bottom panel of Figure 3, and were used for subsequent 
Kepler- related lag analyses. 

We note that the Kepler passband is much better 
matched to the R band rather than the F band, which 
means there could be color-dependent variability signals 
contributing to discrepancies between the F and Kepler 
light curves. We also note that the SAP light curves 
from the Kepler archive are susceptible to several instru- 
mental effects. First, the use of different sized apertures 
between individual quarters affects the SAP fluxes more 
so than the Kepler Pre-search Data Conditioning Simple 
Aperture Photometry (PDCSAP) fluxes because of the 
much smaller aperture sizes of SAP. Additionally, the 
effects of DVA are also larger in the SAP light curves 
compared to the PDCSAP light curves. A more robust 
analysis of Kepler AGN data would require re-extracting 
the SAP light curve over a larger set of pixels to remove 
these systematics, but that is beyond the scope of this 
work. 


3.3. Continuum, Light Curve Characteristics 

To quantify the observed KA1858+4850 continuum 
variability during our monitoring period, we computed 
the statistics ff max and F v ar for consistency with previous 
reverberation mapping studies (Rodriguez-Pascual et al. 
1997; Peterson et al. 2004). ff max is defined as the ratio 
between the maximum and minimum observed fluxes, 
and ff var is defined as 


V a 2 - (s 2 ) 
(f) 


(3) 


where a 2 is the sample variance, (S 2 ) is the mean square 
value of the measurement uncertainties, and (/) is the 
unweighted mean flux. ff var is essentially an estimate of 
the intrinsic root-mean-square (rms) variability relative 
to the mean flux corrected for random errors. We found 
ff max = 1.56 and F var = 0.086 for the F-band light curve 
and -ff ma x = 1-45 and ff var = 0.080 for the Kepler light 
curve. 

A previous AGN monitoring campaign carried out by 
the LAMP 2008 collaboration observed 13 AGNs over 
a two-month period (Bentz et al. 2009). Five of these 
AGNs (Mrk 142, Mrk 1310, Mrk 202, NGC 4253, and 
NGC 4748) are NLS1 galaxies with full width at half- 
maximum intensity FWHM(H/3broad) < 2000 km s -1 . 
The Avar values for their F -band light curves range from 
0.27 to 0.73, and the ff max values range from 1.12 to 
1.39. Compared to these NLSls, KA1858+4850 was sig- 
nificantly more variable during our monitoring period, 
with both ff var and ff max values much larger than those 
for the LAMP 2008 NLS1 galaxies over their two-month 
monitoring period. 


4. SPECTROSCOPIC OBSERVATIONS 

Spectroscopic observations of KA1858+4850 were car- 
ried out using the Kast Spectrograph on the Shane 3-m 
telescope at Lick Observatory. This spectrograph is usu- 
ally mounted only during dark runs. We employed an 
interrupt-mode observing method, where every group of 
Kast observers took one exposure of KA1858+4850 on 
each of their regularly scheduled observing nights. This 
enabled us to spectroscopically monitor the AGN for a 
total of nine months, much longer than what is achiev- 
able by most dedicated observing campaigns at classi- 
cally scheduled facilities. 

The Kast spectrograph has a D55 dichroic that splits 
light from the slit at about 5500 A into separate blue- and 
red-side cameras. Our standard setup used a 600/4310 
grism on the blue side, which gives a wavelength disper- 
sion of 1.02 A pixel -1 and wavelength range of 2090 A. 
However, the wavelength coverage was inconsistent be- 
cause each group used a slightly different blue-side setup 
that shifted the wavelength coverage, and on the nights 
of 2012 February 16, 2012 March 4, 2012 April 19, and 
2012 May 1, the observers employed a 830/3460 grism. 
We used the wavelength range 4000—5500 A for our anal- 
ysis as this is common to all spectra. This wavelength 
range includes the H/3, Hy, Hi5, [O III], and He II emis- 
sion lines, as well as a portion of the Balmer continuum. 

On the red side, because different observing teams used 
significantly different setups for their primary science tar- 
gets, we were unable to obtain a complete set of spectra 
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Fig. 4. — Mean rest-frame spectrum of KA1858+4850 con- 
structed from all nights with both blue- and red-side Kast obser- 
vations. 

with consistent quality and wavelength coverage for anal- 
ysis of the Ha line. For reference, Figure 4 shows the 
unweighted mean AGN spectrum constructed from all 
nights with both blue- and red-side Kast data. 

From February to November of 2012, weather permit- 
ting, each regularly scheduled group of Kast observers 
took at least one 1200 s exposure of KA1858+4850 at 
the parallactic angle (Filippenko 1982) using a 2" slit, 
along with one 120 s exposure of a flux-standard star with 
the same slit width for calibration. Two consecutive ex- 
posures of KA1858+4850 were taken on two nights and 
three consecutive exposures were taken on four nights. 
The flux standards we used are BD+284211, Feige 34, 
G191B2B, and HZ 44, in decreasing order of frequency. 
Very few spectra were taken in February and March ow- 
ing to poor weather conditions. April and May had sev- 
eral good nights of data, and starting from June until 
the end of the campaign in November, we obtained spec- 
tra during more than two-thirds of the Kast nights each 
month. We obtained spectroscopic data from a total of 
74 nights. 

5. SPECTROSCOPIC REDUCTIONS AND 
MEASUREMENTS 

Spectroscopic data reduction included overscan sub- 
traction, flat fielding, cosmic ray cleaning using the 
L. A. Cosmic routine (van Dokkum 2001), extraction with 
a width of 6"88 (corresponding to a 16-pixel extraction 
window for the blue-side data), wavelength calibration 
employing line-lamp exposures, and flux calibration us- 
ing standard stars. We took unweighted extractions for 
AGN spectra and optimal extractions for standard-star 
spectra. Spectra taken on the same night were aver- 
aged into a single spectrum. We also propagated the ex- 
tracted error spectrum through subsequent calibrations 
and analyses. 

We attempted to perform spectral decomposition us- 
ing methods described by Barth et al. (2013) to isolate 
the broad-line components. However, owing to the pres- 
ence of weak and blended emission lines as well as lim- 
ited spectral coverage, many single-epoch spectra pro- 
duced poorly constrained fit parameters for the contin- 
uum components, He II and Fe II emission, and red- 
dening. We therefore used the traditional approach of 
measuring line fluxes by employing a linear fit to ap- 
proximately subtract the continuum underlying emission 
lines. The decomposed components of the higher-SNR 
mean spectrum are displayed in Figure 5 for reference. 



4600 4800 5000 5200 5400 

Observed Wavelength (A ) 


Fig. 5. — KA1858+4850 mean spectrum (black), the combined 
model fit of all components (red), and individual spectral fit com- 
ponents. The H 7 and [O III] A4363 blend was excluded from the 
fits in order to limit the total number of fit parameters. 


To quantify the flux-measurement uncertainty intro- 
duced by using this linear interpolation approach as op- 
posed to the spectral deomposition method, we also mea- 
sured /( H/3) of a series of H/3-only spectra. Each H/3- 
only spectrum was created by subtracting from the data 
all model fit components except the broad and narrow 
H/3 models. We found that /(H/3) measured from the 
H/3-only spectra are, on average, 1.2% higher than those 
measured by simply interpolating over the continuum in 
the data. 

To calibrate the relative fluxes between indi- 
vidual spectra, we followed steps described by 
van Groningen & Wanders (1992), where the [O III] lines 
are taken to be constant in flux for the duration of the 
campaign. As the [O III] line is emitted by gas in the 
narrow-line region, which is much farther out from the 
black hole than the BLR, the time delay in line response 
to continuum variations is much longer than typical re- 
verberation mapping campaigns. The algorithm applies 
a multiplicative flux scaling factor, a small wavelength 
shift, and a convolution with a Gaussian kernel to a re- 
gion in each individual spectrum that contains a nar- 
row emission line and some surrounding continuum, and 
searches for a combination of these parameters that min- 
imizes the residual between this region in the individual 
spectrum and the same region in a reference spectrum. 
We constructed the reference spectrum from the mean of 
all Kast blue-side spectra taken with the 600/4310 grism, 
and chose the observed wavelength range 5390—5410 A, 
which encompasses the [O III] A5007 emission line, to be 
the comparison region. Spectra taken with the 830/3460 
grism were not used to make the reference spectrum, but 
were calibrated using the same method. The flux scale 
factors range from 0.27 to 4.70. The median wavelength 
shift is 1.2 A, which is consistent with the amount ex- 
pected from miscentering the AGN in the slit. 

We followed steps described by Barth et al. (2011) to 
assess the accuracy of the spectral scaling, and calculated 
the normalized excess variance of the [O III] emission- 
line light curve. The normalized excess variance, cr^ x , is 
defined by normalizing Eq. 1 by a factor of the mean flux 
squared, giving 
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TABLE 3 

Wavelength Windows for Flux Measurements 


Line (A) 

Line Window (A) 

Continuum Windows (A) 

H/3 

5200-5290 

5130-5160, 5360-5390 

He II 

4990-5100 

4960-4980, 5120-5160 

h 7 

4650-4720 

4600-4640, 4730-4780 

H<5 

4395-4455 

4360-4380, 4470-4500 

H/3- blue 

5200-5238 

5130-5160, 5360-5390 

H/3-core 

5239-5249 

5130-5160, 5360-5390 

H/3-red 

5250-5290 

5130-5160, 5360-5390 


Note. — Wavelengths are in the observed frame. 


1 N 

-nx=^£[(^-M) 2 -^]. (4) 

We found cr nx « 0.02 for the [O III] light curve after 
flux scaling, indicating that, above the uncertainties from 
photon counting in flux measurements, there is an addi- 
tional scatter on the order of 2% of the total [O III] flux 
in the scaled light curve. This scatter may be caused 
by a combination of variations in seeing, miscentering of 
the AGN in the slit, and nightly variations in the instru- 
ment focus. Overall, this is a relatively small effect on 
the flux scaling of the H/3 light curve. We added this 
2% flux scatter in quadrature to all spectroscopic flux 
uncertainties before performing further analysis. 

The spectroscopic data were photometrically cali- 
brated by carrying out synthetic F-band photometry on 
the spectrum from 2012 September 9, which was taken 
under nearly photometric conditions. We compared this 
magnitude to the aperture photometry magnitude from 
the same night and calculated a scale factor of 1.15 that 
needed to be applied to the spectrum to bring the syn- 
thetic photometry measurement into agreement with the 
aperture photometry measurement. We then applied this 
scale factor to the entire set of Kast spectra. 

To obtain emission-line fluxes, we first subtracted a 
local linear continuum surrounding the line, then inte- 
grated over the emission-line profile. Table 3 shows the 
wavelength ranges used for each line and their local con- 
tinuum windows. Table 4 gives the spectroscopic mea- 
surements of the H/3, Hy, H<5, and He II emission lines for 
the entire dataset, as well as the SNR for each epoch mea- 
sured using the observed wavelength range 4500—4600 A. 
The median SNR per pixel is 28. 

Figure 6 displays the F-band photometric and spec- 
troscopic light curves for the H/3, Hy, H(3, and He II 
emission lines. The scaling routine works best for wave- 
length ranges closest to the [O III] emission lines, so 
at wavelengths farther away from [O III], the higher- 
order Balmer-line light curves become progressively nois- 
ier. Noise in the He II light curve is primarily caused by 
weak line strength as well as a lack of true continuum 
surrounding the line. The presence of Fe II lines blended 
into the blue side of He II, and the fact that the He II line 
is intrinsically very weak and broad, make fitting the true 
continuum with a linear model very difficult. The spec- 
tral decomposition components of He II are also poorly 
constrained owing to the line’s low amplitude. 

Figure 6 also illustrates the spectroscopic light curve 
for the observed wavelength range 4500—4600 A. This 


region is dominated by continuum emission, so its light 
curve can be compared with the F-band light curve. 
This spectroscopic 13-band continuum light curve, de- 
noted by B s , is noisier than that of the F band ow- 
ing to higher susceptibility to seeing variations and slit 
losses, but the two light curves show consistent variabil- 
ity trends during the monitoring period. 

The F var and R m ax values for each of the light curves 
are listed in Table 5. The higher-order Balmer lines 
exhibit distinctly larger relative variability amplitude, 
and the He II line is proportionally more variable than 
all the Balmer lines. Both results are in agreement 
with findings of previous reverberation mapping pro- 
grams (Peterson & Ferland 1986; Dietrich et al. 1993; 
Kollatschny 2003; Bentz et al. 2010). 

Figure 7 shows the mean and rms spectra of 
KA1858+4850 constructed from all blue-side spectra 
taken with the 600/4310 grism after applying [O III] 
spectral scaling. The rms spectrum indicates the amount 
of relative variability at each wavelength. The [O III] 
narrow lines have low residuals in the rms spectrum, in- 
dicating good spectral flux calibration results using the 
[O III] lines. The broad Balmer lines clearly stand out 
with very high variability. He II appears to be highly 
variable in the rms spectrum, even though the line is 
weak in the mean spectrum owing to blending with Fe II 
lines. 

6. LAG MEASUREMENTS 
6.1. Cross-Correlation Measurements 

We calculated the lag between the continuum and 
each emission-line light curve illustrated in Figure 6, as 
well as between the photometric and spectroscopic light 
curves, by employing the interpolation cross-correlation 
technique developed by Gaskcll & Peterson (1987) and 
described by White & Peterson (1994), Peterson et al. 
(2004), and Bentz et al. (2009). We computed the cross- 
correlation function (CCF) for r values from —20 to 40 
days in increments of 0.25 days. The lag for each emission 
line is then calculated in two ways: by using the peak of 
the CCF, defined as Tp ea k, and by using the centroid of 
CCF values above 80% of the peak value, defined as T cen - 
We opted to use r cen for Mbh estimates as Peterson et al. 
(2004) showed that this yields more consistent black hole 
mass estimates between different emission lines. 

In cases where the continuum light curve exhibits dis- 
tinct global trends, a detrending procedure is sometimes 
applied prior to cross-correlation analysis, where a lin- 
ear function is fitted to and subtracted from the light 
curve so that only local variations are taken into account 
in the cross-correlation. We computed the H/3 lag both 
with and without detrending using a linear fit. In the 
case without detrending, the lag uncertainties are smaller 
and the CCF peak is higher, indicating a more robust 
CCF. Therefore, we chose to omit the detrending pro- 
cedure for our final cross-correlation analysis. The top 
panel in Figure 8 shows the CCF for the four emission- 
line light curves with the photometric light curve. We 
also computed the auto-correlation function (ACF) for 
the photometric light curve, which peaks at zero lag as 
expected. 

To determine the final lags and their uncertain- 
ties, we employed the same Monte Carlo bootstrapping 


TABLE 4 

Spectroscopic Measurements for KA1858+4850 


UT Date 

HJD— 2450000 

SNR 

/(H/3) 

/(H 7 ) 

(10“ 15 

/(H<5) 

erg cm 2 s 1 ) 

/(He II) 

2012-02-16 

5974.087 

11 

42.87 ± 0.38 

19.12 ± 0.51 

8.27 ± 0.54 

11.60 ± 0.47 

2012-03-04 

5991.073 

35 

40.90 ± 0.14 

18.37 ± 0.14 

9.84 ± 0.14 

7.30 ± 0.14 

2012-04-02 

6020.024 

19 

42.42 ± 0.30 

19.30 ± 0.27 

12.24 ± 0.27 

9.68 ± 0.28 

2012-04-16 

6033.929 

21 

45.20 ± 0.28 

22.75 ± 0.26 

13.73 ± 0.26 

10.20 ± 0.27 

2012-04-16 

6034.975 

34 

46.41 ± 0.19 

22.53 ± 0.16 

13.15 ± 0.16 

8.04 ± 0.16 


Note. — Listed SNR is the signal-to- noise ratio per pixel for the observed wavelength range 4500-4600 A in the AGN spectra. Measured 
fluxes include the blended broad and narrow emission lines. (Full table available in online version.) 


TABLE 5 

Light Curve Statistics 


Light Curve 

FV ar 

-Rmax 

V 

0.084 

1.50 ± 0.05 

B a 

0.119 

1.88 ± 0.07 

H/3 

0.076 

1.41 ± 0.05 

h 7 

0.078 

1.52 ± 0.07 

H<5 

0.111 

2.19 ± 0.13 

He II 

0.245 

3.54 ± 0.26 


Note. — R ma x and F va r values for V, B s , and the four emission 
lines. Higher-ionization lines show larger variations. 


method used by Barth et al. (2011) and described by 
White & Peterson (1994) and Peterson et al. (2004). We 
constructed 10 4 modified realizations of the continuum 
and emission-line light curves. Each realization is made 
by randomly choosing n data points from the actual light 
curve allowing resampling, where n is the total number 
of points in the dataset. If a point is picked m times, 
then its uncertainty is reduced by a factor of m 1 / 2 . The 
simulated light curves are then varied by adding random 
Gaussian noise based on the measured uncertainties at 
each data point. We then computed the CCF for each 
pair of simulated continuum and line light curves to con- 
struct distributions of r cen values. The median values 
are chosen as the final lag results, and the uncertainties 
on Teen are the la thresholds in the distribution centered 
around the median. 

Table 6 gives the measured r pe ak and r ce n values for 
the four emission- line light curves with respect to the V - 
band light curve. The He II lag is consistent with zero 
within la uncertainties. The larger fractional uncertain- 
ties on the higher-order Balmer line lags, as well as on the 
He II lag, can be attributed to their noisier light curves 
due to less precise spectral scaling at wavelengths farther 
from [O III]. 

The lag times are progressively shorter for higher- 
order Balmer lines. Specifically, we find lag ratios of 
t(H/3):t(H7):t(H 5)=1.00:0.75:0.44. This is consistent 
with the picture of a BLR stratified in optical depth 
(Rees et al. 1989; Korista & Goad 2004), as well as with 
findings from previous reverberation mapping campaigns 
(e.g., Bentz et al. 2010). 

Additionally, we attempted to obtain velocity-resolved 
lag measurements for KA1858+4850, since the lag behav- 
ior as a function of velocity across broad emission fines 
can contain information about BLR kinematics. We di- 
vided the H/3 line profile into three wavelength segments: 
5200—5238 A for the blue wing, 5239—5249 A for the 
core, and 5250—5290 A for the red wing. The H/3 lag for 


TABLE 6 

Observed-Frame Lag Measurements 


Emission Line 

Tpeak (days) 

Teen (days) 

tjavelib (days) 

H/3 vs. V 

o or + 7.25 

o-^0_1.00 

14 . 581 2 ;)® 

14 . 18 !);)! 

H 7 vs. V 

7-50 !);)) 

10 . 9611 ;?® 

10 . 04 !); 4 ! 

H5 vs. V 

6-50 ±1:11 

6 . 421 2 ;®® 

r o-i+l. 06 

O.Oi _2 03 

He II vs. V 

0 7ty + 0 - 50 
u - '°_0.50 

—0 C S8~*~ 1 ' 20 

U.05_o g5 

- 2 - 86 !);!) 

H/3 blue vs. V 

7 oc + 1-00 

' - ZO - 0.75 

13 . 431 );)) 

13 . 85 !);)! 

H/3 core vs. V 

15.00^-75 

15 . 50 !);°) 

14 73+ 0 - 90 

i^t. /O_o 89 

H/3 red vs. V 

8 - 50 t ?;£ 

12 . 89 l 3 3 ;« 

14 . 25 !);)! 

H (3 vs. B s 

11 . 50 ±®;gg 

14 . 89 ! 4 ;)° 

15 . 67 !);)) 

V vs. B s 

2-25 !);)! 

1 ' 68 !l'39 

1 64 +0 ' 30 

- 0.73 

H (3 vs. Kepler 

o oc+6.50 
°‘^- 1.00 

14 17+ 2 - 26 

13 42+ 1 10 

lo.ftz — i IQ 

Hy vs. Kepler 

6-75 ± i;gg 

Q 49+3.02 
v.<±v_2 24 

9- to!);!! 

H<5 vs. Kepler 

4 50"*” 1 ' 75 

^• ou — 1.75 

4 86 +2 ' 78 

ou - 2.27 

4 . 86 !);)! 

He II vs. Kepler 

0 00" 1 ” 0-50 
U.UU—0 75 

—0 72+ 0 - 72 
u - ' 0.72 

U. 55_o 03 

Kepler vs. V 

o-sotS:” 

1 OO -1-0 ' 47 
i.uu_ 0 47 

0-76 !);!) 

Kepler vs. B s 

1 . 751 )®° 

1 95+ 1 - 28 
i.yo_l 16 

2 . 06 !);)! 


Note. — Cross-correlation r pea k, cross-correlation r cen , and 
JAVELIN lags. Observed-frame lags can be converted to rest-frame 
lags by dividing by 1 -j- z. 


each segment is listed in Table 6. We found marginal ev- 
idence for longer lag in the emission-line core and shorter 
lags in the wings. We were unable to obtain useful lag 
measurements for smaller velocity bins, and therefore re- 
frain from drawing any definitive conclusions regarding 
the kinematics of the BLR. 

6.2. JAVELIN 

We used an alternative method of estimating emission- 
line lags, which employs a statistical model for quasar 
variability. This method uses the Python code JAVELIN 
v.0.3a (Zu et al. 2011) to model the optical AGN contin- 
uum variability as a damped random walk process with 
covariance function 

S mw (At) = a 2 exp ^ ^ , (5) 

where r r is the “relaxation time” required for the vari- 
ability to become roughly uncorrelated, and a is the 
variability amplitude on timescales much shorter than 
r r (Kelly et al. 2009). JAVELIN fits r r and a for the 
AGN continuum light curve, then models the emission- 
line light curves as lagged, smoothed, and scaled versions 
of the continuum light curve. An important caveat of 
using JAVELIN for the KA1858+4850 lag analysis is that 
the Kelly et al. (2009) damped random walk model pro- 
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6. — KA1858+4850 V-band magnitude, continuum flux measured from the spectroscopic data, and emission-line light curves. Plotted 
include the 2% flux scatter found by computing the normalized excess variance of the [O III] light curve. 
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Lag (Days) 

Fig. 8 . — Top: Cross-correlation functions between the four emis- 
sion lines and the F-band continuum, and the auto-correlation 
function of the F-band continuum. Bottom: Probability distribu- 
tions of JAVELIN lags for H/3, H 7 , and H<5. These distributions were 
obtained with 6.25 X 10 4 iterations, while the He II distribution was 
obtained with 2.5 X 10 5 iterations. However, the He II distribution 
is poorly constrained owing to the line’s noisy light curve, and is 
therefore omitted in this plot. 

duces variability power spectra with a slope of — 2 , while 
Mushotzky et al. (2011) showed that KA1858+4850 has 
a power-spectrum slope of ~ —3. 

The 14-band light curve was rebinned into one-day in- 
tervals for analysis with JAVELIN in order to cut down 
on computation time. While JAVELIN is, in principle, 
able to fit a large number of emission-line light curves 
simultaneously, the lags were poorly constrained in this 
case for fitting three emission-line light curves simulta- 
neously, most likely because of the monthly gaps in the 
data when the Moon was bright. Therefore, we chose the 



6000 6100 6200 
HJD - 2450000 


Fig. 9.- JAVELIN model results for the continuum (F band), H/3 
(spectroscopic), and H 7 (spectroscopic) light curves (black solid 
lines), the model 1 a uncertainties at each time (shaded regions), 
and observational data and uncertainties. 

two-line analysis method, where we fit each of H 7 , HJ, 
and He II emission-line light curves simultaneously with 
that of H/3. The H/3 lags computed from pairing with 
H 7 and H<5 are consistent with each other, while the H/3 
lag computed from pairing with He II yielded a slightly 
shorter lag. This is likely due to the noisy He II light 
curve as well as the fact that He II intrinsically has a lag 
that is very short compared to the monthly gaps in the 
light curves, which makes the lag difficult to measure. 
We use the H/3 lag value obtained from pairing with H 7 
as tjavelin for H/3. 

Table 6 lists the JAVELIN lags, which are consistent 
with those obtained using cross-correlation techniques 
within ler uncertainties. Lags for the H/3 blue wing, core, 
and red wing were computed simultaneously in a three- 
line JAVELIN run, and the F-band and H/3 lags with re- 
spect to the B s band were obtained from a two-line run. 
The bottom panel of Figure 8 shows the JAVELIN distri- 
butions for H/3, H 7 , and HJ lags, and Figure 9 shows the 
JAVELIN model results for the continuum, H/3, and H 7 
light curves. 

We note that both the CCF and JAVELIN He II lags are 
slightly negative, which is likely caused by the combined 
effects of the higher ionization (and therefore shorter lag) 
of He II, and a slight contaminating lag signal in the V- 
band light curve, described in the next section. 

6.3. H/3 Contamination in the V band 

There is a small biasing factor in the H/3 lag calcula- 
tions from the emission- line contribution to the F-band 
flux. While the F-band light curve would ideally repre- 
sent pure continuum, the presence of the H/3 line in the 
F filter adds a flux component that contains a lag sig- 
nal. Consequently, the calculated lags from the biased 
continuum would be shorter than those obtained with a 
pure continuum. 

To determine the magnitude of this contribution, we 
first combined the blue- and red-side spectra from a sin- 
gle night, and next performed synthetic photometry on 
the spectrum using a Johnson F filter. We then removed 
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the H/3 line from the spectrum by directly interpolating 
over it, performed synthetic photometry on the modified 
spectrum, and compared the two magnitude results. We 
found that H/3 contributes approximately 9.6% of the V- 
band flux, and assume that H/3 dominates the variable 
emission-line contribution in the V band and is therefore 
the main source of the lag bias. 

To quantify the effect of this bias on the calculated 
lag, we simulated 10 4 pure AGN continuum light curves 
using methods described by Timmer & Koenig (1995), 
and simulated corresponding emission-line light curves 
by convolving the continuum light curves with a 6 func- 
tion at a lag of 14 days. We then simulated 10 4 contam- 
inated F-band light curves by adding a lagged emission- 
line contribution to the continuum at the 9.6% level. The 
pure and contaminated continuum light curves in each 
pair are both then cross-correlated with the correspond- 
ing emission-line light curve to create two distributions 
of lag times. We found a median lag of 14.0±2.1 days for 
the pure continuum case and a median lag of 13.2 ± 2.1 
days for the contaminated continuum case, indicating an 
expected bias of 0.8 days. This prediction is similar to 
the bias we find from observations. 

We found an H/3 lag of 14. 89+5 ' 4 q days with respect 
to B s , which contains no H/3 flux contamination, indi- 
cating an observed bias of ~0.3 days compared to the 
lag-contaminated case. We also found a small positive 
lag of 1. 68^1 '39 days for the F light curve with respect 
to B s . However, in both simulations and observations, 
the biases are smaller than the lag uncertainties for H/3 
with respect to the F-band light curve. We therefore 
conclude that the lag bias due to H/3 flux contribution in 
the F band is present, but is small compared to the ler 
uncertainties on the r ce n measurements. 

6.4. Lags with Respect to Kepler Light Curve 

We computed the lag of each emission line with respect 
to the scaled Kepler fluxes using both cross-correlation 
analysis and JAVELIN. The Kepler light curve has a ca- 
dence of 30 minutes, giving a total of ~ 1.3 x 10 4 data 
points over three quarters. We binned the light curve 
into bins of 12 and 24 hours to use in the cross-correlation 
and JAVELIN analyses, respectively. For both CCF and 
JAVELIN, we found the emission-line lags with respect to 
the Kepler light curve, listed in Table 6, to be consis- 
tent with but slightly shorter than those with respect to 
the F-band continuum. This is consistent with expec- 
tations, since the Kepler passband includes Ha, which 
introduces an additional lag signal to the Kepler light 
curve compared to the F-band data. The redder portion 
of the continuum could also have a small lag with re- 
spect to the bluer continuum (Sergeev et al. 2005), since 
the redder continuum emission comes from larger radii 
in the accretion disk than where the F-band continuum 
is emitted. The combined effects of broad emission lines 
and red continuum in the Kepler band should account for 
the shorter emission-line lags measured against the Ke- 
pler light curve as compared to those measured against 
the F-band light curve. 

We also measured the lag of the Kepler light curve with 
respect to both the F-band and J3 s -band light curves, 
and found small positive lags for both cases. This also 
supports the idea of broad emission lines and the red 
continuum introducing a lag signal to the Kepler light 


curve. 

7. LINE WIDTHS AND BLACK HOLE MASS ESTIMATE 

There are two conventional methods of measuring the 
broad-line width: using the FWHM and the line disper- 
sion (u ii n e) of the emission-line profile. The line profile is 
typically taken to be the rms profile, since using the vari- 
able portion of the spectrum instead of the mean spec- 
trum implies a black hole mass estimate based only on 
components of the emission line that echo the contin- 
uum signal (Peterson et al. 2004). The line dispersion is 
defined as 


CT line — 


MVem 
J { E Si 



(6) 


where Si is the flux density at wavelength bin A i and Ao is 
the flux-weighted centroid wavelength of the line profile. 
In this empirical method of measuring the line width, 
the line profile is not fitted to any functional model. We 
used the same line and continuum windows to measure 
the line width as those used in measuring line fluxes. 

To determine the final FWHM and cru ne values and un- 
certainties, we employed the bootstrap method described 
by Peterson et al. (2004). The entire dataset contains N 
spectra. For each bootstrap realization, we randomly se- 
lected N spectra from the dataset allowing reselection, 
constructed the mean and rms line profiles from this ran- 
domly sampled set, and measured the line dispersion of 
the rms profile. From multiple realizations, we built up 
a distribution of FWHM and line-dispersion values, and 
took the median and standard deviation of the distribu- 
tions to be the final FWHM and crime and their uncer- 
tainties, respectively. We removed the instrumental line 
width by taking the width of the A5086 Cd I calibration 
line in a 2"-slit width exposure and subtracting it from 
the measured FWHM or crii ne in quadrature. We found 
FWHM = 324 km s _1 for the Cd I calibration line for a 
Gaussian fit to the line profile. After correcting for the 
instrumental line width, we found FWHM = 1511 ± 68 
km s” 1 and cri ine = 770 ± 49 km s _1 for the H/3 line in 
the rms spectrum. 

We also measured the H/3 FWHM and dime for the 
mean profile. To ensure exclusion of the narrow-line 
component in the width measurements, we measured the 
FWHM and CTii ne of the broad H/3 model based on the 
spectral decomposition of the mean spectrum, as shown 
in Figure 5. The [O III] narrow-line profile was used 
to model the narrow H/3 line in the spectral fitting rou- 
tines, and / (H/3 narrow ) / / ( [O III]a 5007) ~ 0.09. We also 
measured the FWHM and di; ne of the broad H/3 model 
for each epoch in our dataset, and took the standard 
deviations about the means to be the FWHM and cr line 
uncertainties. We found FWHM = 1820 ±79 km s _1 and 
cr line = 853 ± 34 km s” 1 for H/3 in the mean spectrum 
after subtracting the instrumental line width. This is 
consistent with previous findings that line widths mea- 
sured from mean spectra tend to be larger than those 
measured from rms spectra (e.g., Bentz et al. 2009). 

The reverberation lag and line width of H/3 combined 
can give a virial estimate of the central black hole mass, 
given by 
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M B h 


(ct)(AH) 2 

G 


(7) 


where r is the H/3 lag time with respect to the continuum 
and ct gives the mean radius of the BLR, AV is the H/3 
line width, G is the gravitational constant, and / is a 
scaling factor of order unity that depends on the incli- 
nation and kinematics of the BLR. Traditionally, since 
these properties of the BLR are usually unknown, the 
scale factor / is chosen to be a value that brings the set 
of reverberation mapped AGNs into agreement with local 
quiescent galaxies in the Mbh — cr* relation, which relates 
black hole mass to host-galaxy bulge stellar velocity dis- 
persion (Onken et al. 2004; Woo et al. 2010; Park et al. 
2012; Grier et al. 2013). 

We use CT cen for the BLR radius and H/3 line dis- 
persion a ij ne of the rms line profile for AH, for con- 
sistency with Peterson et al. (2004), and a scale factor 
of / = 5.13, calculated by Park et al. (2012) based on 
the updated local AGN Mbh — <7* relation obtained 
with the forward regression method. Combining the 
H/3 vs. V lag of i4.58l2.50 days, corresponding to a 
rest-frame lag of r cen = 13.53/232 days, and cr line = 
770 ± 49 km s -1 , we obtain a virial black hole mass 
estimate of Mbh = 8.O6II72 x 10 6 M 0 . If we follow 
the prescription of Grier et al. (2013) by using the H/3 
Tjavelin and a scale factor / = 4.31, we find M B h, javelin = 
6-58lJ;98 x 10 6 M 0 . 

The above uncertainties on Mbh include only errors 
propagated from the lag and emission-line-width mea- 
surements. If we incorporate the uncertainties on the 
mean scale factor from the linear fits by Park et al. 
(2012), / = 5.13 ± 1.30, then our black hole mass esti- 
mate becomes Mbh = 8.O6I2.6? x 10 6 AL 0 - ^ ' s evident 
that true uncertainties on the virial estimate of Mbh are 
dominated by the systematic uncertainties in the scale 
factor, which are significantly larger than those derived 
from the lag and line-width measurements alone. 

We note that there are other estimates of the scale fac- 
tor, such as those obtained by separating galaxies into 
different populations based on mass (Greene et al. 2010) 
and morphology (Graham et al. 2011), which yield esti- 
mates of / different from that of Park et al. (2012) by 
up to a factor of ~2. For example, Woo et al. (2013) in- 
vestigated the scale factor for both quiescent and active 
galaxies as a combined sample and found / = 5.9j/L 5 . 
Furthermore, recent work by Ho & Kim (2014) showed 
that the scale factor can be different for galaxies with 
pseudobulges and classical buldges, with / = 3.2 ±0.7 for 
pseudobulges and / = 6.3 ± 1.5 for classical bulges. Var- 
ious ongoing efforts that further examine the Mbh — cr* 
relation for local galaxies will improve the precision of 
the scale factor in the near future as the number of re- 
verberation mapped AGNs increases. Moreover, there 
has been progress in constraining / for individual galax- 
ies by dynamically modeling the BLR (Pancoast et al. 
2013). 

In addition, we obtained Mbh estimates using the 
broad H7 and H(5 lines. No lag estimate was attempted 
using He II since the line has a negative lag. The line 
widths, rest-frame lags, and derived Mbh values are 
listed in Table 7. The H7 and Hd light curves are signif- 


icantly noisier than that of H/3; thus, it is not surprising 
that, for both CCF and JAVELIN cases, the derived Mbh 
values have much higher fractional uncertainties com- 
pared to the H/3 Mbh- For both CCF and JAVELIN lags, 
the H7 Mbh estimates, though consistent with the H/3 
Mbh values within la uncertainties, are slightly smaller 
than those of H/3, and Mbh estimates for Hd are smaller 
still. This may be due to the fact that we are using the 
same / factor for all the emission lines, while the strati- 
fied nature of the BLR may imply different scale factors 
for each line that depend on the geometry and kinematics 
of the line-emitting gas. 

We would like to compare KA1858±4850 to other 
AGNs having similar black hole masses by studying its 
location on the Mbh — cr* relation as well as the Mbh — 
Lbuige relation (black hole mass vs. host-galaxy bulge 
luminosity). However, because KA1858±4850 appears 
point-like at ground-based resolution, it is impossible to 
observe structural properties of the host galaxy without 
high-resolution images from the Hubble Space Telescope 
(HST) or ground-based observations using adaptive op- 
tics. Additionally, our Lick spectra cannot be used to 
measure stellar velocity dispersion in KA1858±4850 ow- 
ing to the galaxy’s weak starlight component compared 
to its AGN luminosity. 

8. EDDINGTON RATIO 

AGNs have been observed to follow a tight correlation 
between the size of the BLR (7?blr) and continuum lu- 
minosity ( L \ ). The 7?blr — L\ relation can be written 
in the form 


log 


/ -Rblr \ 
\1 lt-day / 


K + a log 


ALa 


10 44 erg s 


(8) 


where L\ is measured at A res t = 5100 A. Bentz et al. 
(2013) found the values of K and a to be 1.560/to!o24 an d 
0.546lg Q27 respectively, with a scatter of around 0.13 
dex for their best fit. From these parameters and the 
lag for KA1858±4850, we expect to find ALa(5100 A) = 
1.64j/g;g5 x 10 43 erg s -1 . We used combined Kast blue- 
and red-side spectra to measure L\ and adopted the spec- 
tral fitting components for the mean spectrum (shown in 
Figure 5) to estimate the starlight contribution in this 
region, which we found to be approximately 40% of the 
total flux. Correcting for Galactic extinction, we roughly 
estimate A/a (5100 A) ss 1.6 x 10 -12 erg cm -2 s -1 for the 
AGN, corresponding to ALa(5100 A) ss 2.4 x 10 43 erg s -1 
for a luminosity distance of 354 Mpc. (We assume the 
same standard ACDM cosmology as Bentz et al. 2013, 
where Hq = 72 km s -1 Mpc -1 , Ou = 0.3, and Ha = 0.7.) 
This is consistent with expectations given the scatter in 
the fit values from Bentz et al. (2013). 

KA1858±4850 is a NLS1, a class of objects thought to 
have high L/Le<m (Pogge 2011, and references therein). 
We apply the bolometric correction used by Kaspi et al. 
(2000), where Lboi ~ 9ALa(5100 A), and obtain an es- 
timate of Lboi = 2.2 x 10 44 erg s -1 and L/Le<m ~ 0.2 
using Mbh = 8.06i)' 3 2 x 10 6 M 0 . We compared this 
Eddington ratio to those of the four LAMP 2008 NLS1 
galaxies, which were calculated using black hole masses 
published by Bentz et al. (2009) and ALa(5100 A) values 
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TABLE 7 

Line Widths, Lags, and Derived Black Hole Masses 


Emission Line 

(Tline (km S U 

Lag Computation Method 

Teen, rest (days) 

AIbH (10 6 Mq) 

H/3 

770 ± 49 

CCF 

1 Q co+2.03 
io.oo_2 32 

8-06ti;SS 



JAVELIN 

IQ ik + 1-08 

TO. ±<J _ 1 00 

6.85lJ;“° 

H 7 

741 ± 73 

CCF 

10. 17^2.56 

5.59i);? 2 



JAVELIN 

9-3i±ti? 

Q qq+ 0-97 
o.yy_Q 98 

H<5 

827 ± 83 

CCF 

5-96 tUl 

4-2o±i:iS 



JAVELIN 

5.39l°;®* 

3-19t°;! 


Note. — Line lags are measured against the U-band continuum. Mbh from CCF lags were calculated using / = 5.13 (Park et al. 2012), 
and Mbh from JAVELIN lags were calculated using / = 4.31 (Grier et al. 2013). 


given in Bentz et al. (2013). After applying the same 
bolometric correction as for KA1858+4850, we found 
L/L e dd = [0.5, 0.7, 1.2, 0.9] for Mrk 1310, Mrk 202, NGC 
4253, and NGC 4748, respectively. Compared to these 
NLSls, KA1858+4850 has a significantly lower Edding- 
ton ratio. 

A recent study by Du et al. (2014) measured the H/3 
lag and Mbh of three NLS1 galaxies (Mrk 335, Mrk 142, 
and IRAS F12397), all of which appear spectroscopi- 
cally similar to KA1858+4850. The authors compute 
the Eddington rate based on a thin accretion disk model 
(Shakura & Sunyaev 1973). This rate, denoted by ro ss , 
is written as 

m ss « 20.1 ( — — ) Mf 2 r/ ss , (9) 

\cos i ) 

where they define L44 = ALa/10 44 erg s _1 at A = 
5100 A, M 7 = M bh / 10 7 M q , and cos i = 0.75 as the 
inclination typical of Type 1 AGNs. For a minimal ra- 
diative efficiency of 7/ ss = 0.038, they find Eddington ra- 
tios of 0.6, 2.3, and 4.6 for Mrk 335, Mrk 142, and IRAS 
F12397, respectively. Applying this same prescription to 
KA1858+4850, with M 7 = 0.81 and L44 = 0.24, we find 
rh ss = 0.2, in agreement with our L/ Xe<m value obtained 
using the Kaspi et al. (2000) bolometric correction. 

9. SUMMARY 

We photometrically and spectroscopically monitored 
the Kepler - field AGN KA1858+4850 over a period of 
nine months. We found an H/3 rest-frame lag of 
I3.53I2 32 days with respect to continuum variations us- 
ing cross-correlation methods, and a lag of 13. lSjli'.Qo 
days using the JAVELIN method. We also measured 
emission-line lags with respect to the Kepler light curve 
and found slightly shorter lags compared to those mea- 
sured against the E-band light curve, which is expected 
given the contributions of broad emission lines and red 
continuum flux to the Kepler band. We measured an 
H/3 velocity dispersion of CTii ne = 770 ± 49 km s _1 , and 
calculated a black hole virial mass of Mbh = 8.06L' 7 2 x 
10 6 Mq using t ccf and scale factor empirically derived 
from local active galaxies by Park et al. (2012), and a 
black hole mass of Mbh, javelin = 6.58tagg x 10 6 M 0 
using r javelin and scale factor taken from Grier et al. 
(2013). For this mass, the Eddington ratio is L/TecM ~ 
0 . 2 . 

KA1858+4850 was the second AGN for which data 
was obtained in this interrupt observing mode from Lick 


Observatory, and the second AGN in the Kepler field 
to be monitored by ground-based telescopes (the first 
being Zw 229-015, Barth et al. 2011). Comparing our 
lag results with those obtained by the LAMP 2008 col- 
laboration (Bentz et al. 2009), our lag uncertainties are 
slightly larger. However, considering the much longer lag 
of KA1858+4850, our H/3 fractional lag precision, at less 
than 20%, is still very good. Our analysis using Kepler 
light curves also offers one of the first direct comparison 
of reverberation mapping results between ground- and 
space-based observations for a Kepler AGN. The suc- 
cess of our campaign demonstrates the robust capabili- 
ties of interrupt-mode observations for monitoring AGN 
variability. Factors that negatively impact our measure- 
ments, such as inconsistency in data quality and gaps 
in the spectroscopic light curves (in this case due to the 
AGN being observed only during dark runs), are miti- 
gated by a well-sampled E-band light curve obtained by 
combining observations from several ground-based tele- 
scopes as well as the long duration of the program. 

Further observations of KA1858+4850 can provide 
additional insight into various properties of the AGN 
and host galaxy. Specifically, observations of the bulge 
properties can put KA1858+4850 on the Mbh — Lbuige, 
Mbh — Mbuige) and Mbh — cr* relations. Since the host 
galaxy is very compact, high- resolution HST or adaptive 
optics imaging will be needed to examine the host-galaxy 
morphology. 
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